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We  present  an  ab  initio  density-functional  calculation  of  the  interaction  of  hydrogen  with  0d 
the  Pd(001)  and  Pd(110)  surfaces  and  with  bulk. palladium.  Our  results  for  the  H-Pd  system  :tlon_ 

are  complemented  by  calculations  for  the  Hj  molecule,  the  bulk,  and  the  clean  (001)  and  (110) - - 

surfaces  of  Pd.  In  the  case  of  H  in  bulk  Pd,  we  determine  the  equilibrium  geometry  and 
electronic  structure  and  discuss  the  nature  of  the  Pd — H  bond.  For  the  adsorption  system 
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H/Pdf001  land  H/?d(110),  we  determine  the  preferential  adsorption  sites,  bond  lengths,  and 
vibration  frequencies  of  the  adsorbate.  We  relate  H-induced  softening  of  Pd-surface-phonon  llty  Cod03 
frequencies  to  a  corresponding  softening  of  Pd-Pd  interactions  in  the  bulk  hydride  phase.  .1  and/or 
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I.  INTRODUCTION 

The  interaction  of  hydrogen  with  transition  metals  is 
both  fundamentally  and  technologically  important.1  The 
dissociative  adsorption  of  molecular  hydrogen  at  a  metal 
surface  and  the  subsequent  surface  and  bulk  diffusion  of 
atomic  hydrogen  is  a  prototypical  surface  process.  This  is 
a  well-suited  model  system  for  the  study  of  the  hydrogen- 
metal  bond  at  different  stages  of  the  reaction.  At  least 
of  equal  importance  is  the  microscopic  understanding  of 
the  effect  of  hydrogen  on  the  electronic  and  structural 
properties  of  the  substrate  metal.  This  effect,  which  has 
received  less  attention  in  the  literature,  manifests  itself 
in  H-induced  surface  relaxation  and  reconstruction,  in 
phonon  softening  and  decohesion,  as  well  as  significant 
changes  of  the  electronic  structure  occurring  both  at  the 
surface  and  in  the  bulk  of  the  substrate  metal. 

The  dissociation  of  Ho  at  transition-metal  surfaces  is 
an  important  catalytic  reaction  that  plays  a  key  role  in 
olefin  hydrogenation,  the  Fischer-Tropsch  Drocess,  and 
the  ammonia  synthesis.2  The  presence  of  hydrogen  in 
bulk  metals  causes  material  embrittlement  and  cracking.2 
On  the  other  hand,  bulk  metals  are  ideal  media  for 
safe  storage  of  hydrogen  in  large  quantities.  More  re¬ 
cently,  substantial  scientific  effort  has  been  made1  5  to 
disprove  reports  of  an  apparent  "cold-fusion'  reaction0 
in  deuterium-loaded  transition  metals. 

We  select  the  hydrogen-palladium  system  for  our 
study,  since  it  exhibits  a  -ch  variety  of  the  above- 
fti  descriDed  pneuomena  that  closes  this  system  to  be  of 
substantial  technological  interest.  The  physical  proper¬ 


ties  of  hydrogen  on  the  Pd(00i)  surface  have  been  inves¬ 
tigated  extensively  by  a  variety  of  experimental  meth¬ 
ods  such  as  thermal  desorption  spectroscopy  (TDS),,  d 
work-function  (A0)  measurements,  and  low-energy  elec¬ 
tron  diffraction  (LEED),8  further  by  He  scattering,9 
electron  energy-loss  spectroscopy,1011  and  transmis¬ 
sion  channeling.12  Hydrogen  adsorption  on  this  sur¬ 
face  has  been  investigated  previously  bv  first-principles 
calculations13  for  different  coverages.  Equilibrium  ad¬ 
sorption  of  hydrogen  on  the  (110)  surface  of  Pd  has 
been  investigated  by  LEED,14  low  energy  ion  scatter¬ 
ing  (LEIS),15  He  diffraction,16  high-resolution  electron- 
energy-loss  spectroscopy  (HREELS), 17,18  work-function 
and  thermal  desorption  spectroscopy,19,20  and  Ly  nu¬ 
clear  reaction  analysis  and  Rutherford  backscattering.20 
LEED  studies  indicate  that  at  hydrogen  coverages  0  > 
1.5,  the  Pdf  110)  surface  exhibits  a  (1  x  2)  reconstruction, 
probably  of  the  missing-row-type,21-23  even  though  the 
atomic  positions  in  the  substrate  remain  somewhat  con¬ 
troversial. 

The  limited  amount  of  theoretical  information  avail¬ 
able  is  caused  by  the  intrinsic  difficulty  to  treat  extended 
transition-metal  surfaces  owing  to  the  coexistence  of  lo¬ 
calized  d  electrons  and  delocalized  sp  electrons.  Use¬ 
ful  information  about  the  Pd-11  interaction  can  be  ob¬ 
tained  from  accurate  quantum-chemical  calculations  for 
very  small  Pd„IIm  clusters24  or  parametrized  calcula¬ 
tions  for  extended  surfaces  based  on  the  embeaoeu-atom 
method,25,26  the  embedded-cluster  method2',  the  effec¬ 
tive  medium  method,28  or  the  tight-binding  formalism.-9 
Only  recently,  however,  have  suificiently  powertul  meth- 
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ods  been  introduced30- 34  that  describe  adsorbates  on 
extended  surfaces  with  adequate  accuracy.  We  report 
here  an  ab  initio  calculation  of  the  hydrogen-palladium 
interaction  at  the  surface  and  in  the  bulk.  Some  of  our 
results  for  H/Pd(00i)  and  II  in  bulk  Pd  have  been  re¬ 
ported  eisewher.13  5  Our  calculations  for  the  interaction 
of  H  with  the  (001)  and  (110)  surface  and  the  bulk  of  Pd 
include  equilibrium  geometries,  bond  energies,  vibration 
frequencies,  and  the  electronic  structure  and  the  nature 
r  of  the  H-Pd  bond.  These  studies  are  completed  by  cor¬ 
responding  calculations  for  the  H?  molecule,  as  well  as 
_ hydrogen-free  Pd(bulk),  Pd(001),  and  Pd(110). 

The  paper- is  structured  as  follows.  In  Sec.  II.  we  out¬ 
line  the  application  of  the  ab  initio  density-functional  for¬ 
malism  to  the  H/Pd  system.  In  Sec.  Ill,  we  present  our 
results  for  the  pure  constituents,  namely  hydrogen,  bulk 
Pd,  bulk  PdH,  and  PdH?  as  well  as  H-free  Pd(001)  and 
Pd(110)  surfaces.  In  Secs.  IV  and  V,  we  describe  the  H- 
Pd  interaction  with  Pd(001)  and  Pd(110),  respectively. 
Finally,  in  Sec.  VI,  we  summarize  our  results  and  present 
conclusions. 

II.  THEORY 

We  use  the  local-density  approximation  (LDA)  (Ref. 
35)  to  determine  the  total  energy,  the  charge  distribution, 
and  the  band  structure  of  the  Pd-H  system  in  the  ground 
state.  The  electron  density  p(r)  is  obtained  by  solving 
self-consistently  the  Kohn-Sham  equations35 

(-7—  V3  +^00  +  V'h  +  (!) 

_  and 

occ 

P(r)  =  2Zl^k(r)l2'  (2) 

-  nk 


into  one  entity. 

_ Then,  the  total  energy  is  given  by 

occ 

£\ot  =  t nit  -  ^  / dr  vh{o( r))p(r) 

-  nk 

-  J  dr  Vxc(p(r))p(r)dr 

4-  J  dr  <xc{p(r)}o'  ’•)  —  .  (3) 

-  where  f.„c{p(r))p\  v)  is  the  exchange-correlation  energy 

density  anti  E,on_,on  describes  the  interaction  between 
. —  1  lon  cores.  In  LDA,  !,c[p(i-)]  is  approxnnateu  by  a  local 
function  t  h;em(pln).35 

_  The  use  ol  local  orbitals  in  the  pseudopotontiai  local 


orbital  method34  makes  it  computationally  feasible  to  use 
high-energy  cutoffs  in  the  seif-consistent  calculation  of 
Fourier  components  of  the  charge  density  and  the  poten¬ 
tial.  Our  basis  consists_of  local  Gaussian-cvpe  orbitais  of 
the  form 

faim(  r)  =  .da|me  a  A|m(#,<J>),  (4) 

which  are  localized  on  atomic  sites  and  can  be 
augmented  by  floating  orbitals.  Adm  are  normal¬ 
ization  constants  and  A'im  axe  the  Kubic  harmon¬ 
ics  {1  ,z,y,z,x2,y2,-2  ,zy,yz.zx}  for  1=0,1, 2.  We  use 
norm-conserving  ionic  pseudopotentials  of  the  Hamann- 
Scliliiter-Chiang-tvpe36  and  the  Hedin-Lundqvist37  form 
of  the  exchange-correlation  potential  in  the  local-density 
approximation. 

Our  criterion  for  the  self-consistent  solution  of  Kohn- 
Sham  equations  is  a  requirement  that  the  maximum 
change  in  any  Fourier  component  of  the  total  potential 
<5V(G)  between  two  successive  iterations  is  smaller  than 
10-4  Ry.  The  Brillouin  zone  is  sampled  by  a  homoge¬ 
neous  mesh  of  k  points  that  are  generated  using  a  special- 
point  scheme.38  Within  the  limits  of  the  computational 
technique  and  the  finite  basis  set,  the  typical  precision  of 
the  total  energy  is  of  the  order  <10-4  eV. 


III.  PURE  CONSTITUENTS 
A.  Hydrogen 

Pseudopotentials  for  isolated  hydrogen  atoms  or 
molecules  are  introduced  for  the  sake  of  convenience 
only.  If  the  cutolT  radius  for  the  pseudopotential  is  small 
enough  (in  our  calculation  rc  =  0.35  a.u.),  the  eigenvalues 
and  eigenstates  are  insensitive  to  the  singularity  in  the 
potential  at  the  origin.  Using  the  Hedin-Lundqvist  form 
of  the  exchange-correlation  potential,  we  find  the  total 
energy  of  the  hydrogen  pseudoatom  to  be  -13.26  eV, 
talcing  into  account  the  calculated  spin-polarization  en¬ 
ergy  of  1.04  eV. 

For  hydrogen  in  a  crystal  environment,  we  use  a  Gaus- 
,  sian  basis  consisting  of  s  and  p  states  with  four  radial 
decays  each,  i.e.,  16  independent  functions.  The  decay 
constants  a,  obtained  by  mimtViizing  the  total  energy  of 
hydrogen  and  fitting  the  ground-state  wave  function,  are 
0.20,  0.42,  0.68,  and  1.36  in  atomic  units.  For  the  calcula¬ 
tion  of  H?,  an  arrangement  of  superceils  with  the  dimen¬ 
sions  9  a.u.xO  a.u. x  11  a.u.  has  been  chosen  to  minimize 
wave-function  ovei.ap  oetween  neighboring  cells.  We  use 
a  cutoff  energy  of  64  Ry  for  the  Fourier  components  of  the 
potential  and  charge  density.  In  the  'E*  ground  state 
z:  th :  ,T .  Jte,  wo  obtain  .%  —  3.7S  A  for  the  equi¬ 

librium  distance  (expt.39  0.74  A.)  and  u,  =  4100  cm-1 
for  the  vibration  frequency  (expt.39  4  4  01  cm-1).  A  dis¬ 
sociation  energy10  Df  =  4.71  eV  (expt.39  D .  =  4.75  eV) 
is  obtained  by  comparing  the  total  energy  of  II?  to  that 
of  a  superceil  arrangement  of  hvdrogen  atoms  (corrected 
liv  the  spm-polarizauon  energy).41 


Here,  V\on  is  the  potential  due  to  the  ion  cores,  Vn  is 
the  Hartree  potential,  and  Vxc  is  the  exchange-correlation 
potential  for  the  charge  density  p(r).  In  the  pseudopo- 
tential  formalism,  we  consider  only  the  valence  charge 
density,  and  replace  V,on  by  a  first-principles  pseudopo¬ 
tential  that  combines  the  nucleus  and  the  core  electrons 


2. 


B.  Palladium 


Since  electronic  and  structural  properties  of  Pd  have 

—  been  studied  previously, 42  we  quote  only  those  results 
that  are  relevant  for  a  comparison  with  hydrogen  com¬ 
pounds.  We  find  the  energy  of  the  isolated  Pd  pseu- 

—  doatom  to  be  796.414  eV,  using  the  lledin-Lundqvist, 
form  of  the  exchange-correlation  potential.  We  consider 
s,  p,  and  d  valence  orbitals  with  four  radial  Gaussian 

—  decays  a  each,  i.e.,  40  independent  basis  functions  for 
each  site.  The  palladium  decays,  which  minimize  the 
bulk  total  energy,  are  a  =  0.27,  0.77,  2.18,  and  6.20  in 

—  atomic  units.  Equilibrium  properties  of  bulk  fee  Pd  are 
obtained  by  considering  Fourier  components  of  the  po¬ 
tential  and  charge  density  up  to  a  cutoff  energy  of  64  Ry. 

—  We  sample  the  irreducible  part  of  t.he  Hulk  Brillouin  zone 
by  47  k  points.  Our  resu  _>  for  the  calculated  cohesive  en-  : 
ergy  £coh  of  bulk  Pd  as  a  function  of  the  lattice  constant  : 

—  a  are  shown  in  Fig.  1.  The  calculated  equilibrium  lac-  •’ 
Cice  constant  a  =  3.90  A  is  in  excellent  agreement  with  i 
the  experimental  value'13  of  3.89  A.  The  bulk  modulus  : 

-  of  2.15xl012  dyn/cm2  agrees  well  with  the  experimen-  i 

tal  value43  of  1.81  x  1012  dyn/cm2.  The  calculated  to¬ 
tal  energy  at  equilibrium  of  —799.79  eV  corresponding  j 

- to  a  Pd(buik)  cohesive  energy  of  3.37  eV  (with  respect  i 

to  the  isolated  pseudoatom)  underestimates  the  expert-  • 
mental  value43  of  3.89  eV.  A  possible  reason  for  the  low 

-  value  of  the  Pd  cohesive  energy  is  an  insufficient  basis  j 

set.  Actually,  we  do  find  that  extending  our  basis  set  by  i 
floating  orbitals  at  interstitial  sites  increases  the  vaiue  of  t 

—  the  calculated  cohesive  energy  by  0.2  eV,  but  does  not : 
affect  energy  differences  that  occur  in  calculations  of  elas-  •' 
tic  properties,  etc.  We  also  checked  the  convergency  of  1 

-  our  results  with  respect  to  the  cutoff  energy.  Due  to  the  i 

smooth  Pd  pseudopotential  that  we  use,  Fourier  compo-  ; 
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—  a  (A) 

—  FIG.  1.  Cohesive  energy  £<■  ,n  .is  a  function  of  the  lattice 
constant  a  for  bulk  Pd  (solid  line)  and  bulk  1**1  II  (dashed 

_ line). 


FIG.  1.  Calculated  electronic  density  of  states  :V(£)  of 
(a)  Pd.  (b)  PdH.  and  (c)  PdH?  (JH-h  =  0.74  A). 


FIG.  3.  ta(  Conventional  unit  cell  of  Pu(fcc),  with  the 
octahedron  cage  shown  by  a  solid  line.  Calculated  charge 
density  p( r)  of  (b)  bulk  Pd.  (c)  Pdil,  and  (d)  Pdllj  (du-ii 
=  0.74  A).  All  contour  plots  are  shown  in  the  (1001  plane, 
which  contains  the  four  circled  Pd  atoms  of  the  octahedron 
cage.  The  position  of  Pd  atoms  is  marked  bv  •:  that  of  II 
atoms  bv  x.  The  charge  density  is  m  units  of  U)"£^ h.u.3  ;  the 
ratio  of  two  consecutive  contours  is  u  -*•  I )//;(  n  )=  1 .5. 


nents  of  the  charge  density  die  off  rapidly  and  the  total 
energies  obtained  with  a  64-  and  190-R.v  cutoff  differ  by 
less  than  0.1  eV. 

The  electronic  structure  of  Pd  is  characterized  by  a 
nearly  filled  4d  band  and  a  partially  filled  5s  band.  The 
density  of  states  N(E)  of  Pd.  presented  in  Fig.  2(a), 
shows  the  Fermi  level  near  the  top  of  the  5-eV-wide  4d 
band.  The  large  value  of  .V(Ef)  is  responsible  for  very 
good  conduction  properties  of  this  metal.  A  schematic 
drawing  of  the  Pd  structure,  depicting  the  octahedral 
cage,  is  given  in  Fig.  3(a).  The  charge  density  of  Pd  in 
the  (001)  plane  (containing  four  octahedral  cage  atoms) 
—  is  shown  in  Fig.  3(b). 

_  C.  Palladium  hydride  and  palladium  dihydride 

Palladium  hydride  PdH  is  a  stable  compound  that 
forms  spontaneously  by  exposing  palladium  metal  to  a 
hydrogen  atmosphere.  In  this  process,  the  palladium  lat¬ 
tice  retains  the  fee  structure,  but  the  lattice  constant  is 
observed  to  increase  by  6%  due  to  hydrogen  uptake  in 
the  high  density  phase.'14  H  atoms  occupy  the  octahedral 
interstitial  sites  in  the  Pd  lattice  [see  Fig.  3(a)],  which 
corresponds  to  the  NaCl  structure  for  the  compound. 

In  an  attempt  to  understand  the  equilibrium  prop¬ 
erties  of  Pd  at  various  concentrations  of  hydrogen,  we 
present  results  for  both  PdH  and  a  hypothetical  com¬ 
pound  PdH2-  The  latter  results,  interpreted  for  the 
chemically  equivalent  compound  PdD;>,3  have  implica¬ 
tions  on  the  feasibility  of  the  reported  "cold-fusion”  re¬ 
action  in  a  Pd  matrix.0 

Basis  functions  for  Pd  and  II  have  been  discussed  in 
the  preceding  subsections  on  the  elements.  In  analogy 
to  bulk  Pd,  we  use  an  energy  cutolf  of  64  Ry  in  the 
Fourier  expansion  of  the  charge  density  and  potentials  in 
order  to  ensure  complete  convergence  of  the  LDA  spec¬ 
trum  and  total  energies.  The  LDA  charge  density  and 
potentials  have  been  obtained  by  sampling  the  Brillouin 
zone  with  a  fine  mesh  of  1331  k  points.  Our  total-energy 
calculation  for  PdH  yields  the  equilibrium  lattice  con- 
stant  a  =  4.1  A,  which  compares  well  with  the  experi- 

_ mental  value44  4.12  A.  The  calculated  bulk  modulus  of 

1.95xl012  dyn/enrr  lies  =sl0%  below  the  value  of  bulk 
Pd,  indicating  H-induced  softening 
_____  We  define  the  cohesive  energy  of  PdIIn  with  respect  to 
isolated  atoms  by 

£coh(P<iHn)  =  £’toiBi(PdH„)  -  EtoiadPd  atom) 

~~  —  n£totai(H  atom)  .  (5) 

For  PdH,  we  find  /TcohlPdll)  =  -6.2S  eV,  which 
- compares  favorably  with  the  experimental  value45 

(PdH)  =  -6.65  eV. 

In  order  to  obtain  the  equilibrium  structure  of  PdlU, 
_  we  first  study  the  interaction  potential  between  two  hy¬ 
drogen  atoms  sharing  the  same  octahedron  cage.5  A  de- 

_ crease  of  the  hydrogen  (or  deuterium)  bond  length  due 

to  the  metal  matrix  is  likely  to  favor  a  "cold-fusion”  re¬ 
action  between  deiitcrons.  Such  a  bond-length  reduction 


is  conceivable  in  the  case  that  II?  aligns  with  its  axis 
between  the  bottom  and  the  top  Pd  aLoms  in  the  octa¬ 
hedron  cage.  Our  results  indicate  that  in  this  geometry, 
the  intramolecular  distance  d(H?)  is  expanded  to  0.94  A 
in  the  Pd  lattice.  This  makes  the  "cold”  nuclear  fusion  in 
Pd  even  less  probable  than  in  the  D?  gas  phase.  For  this 
equilibrium  H?  bond  length,  we  find  the  cohesive  energy 
of  PdH?  to  be  £’con(PdHo)  =  -6.7S  eV. 

Our  results  have  been  confirmed  by  independent  cal¬ 
culations  of  Ho/Pdfbulk),  which  also  consider  different 
molecular  orientations.46  These  calculations  predict  a 
spontaneous  reorientation  of  the  Ho  molecule  from  the 
[001]  to  the  [111]  direction,  accompanied  by  a  further 
increase  of  the  bond  length,  y  . , 

The  stability  of  palladium/hydnde  compounds  can  be  ^u'1 
investigated  best  by  comparing  formation  energies  of 
PdH„  (starting  from  bulk  palladium  and  molecular  hy¬ 
drogen),  which  is  given  by 


~form(PdH„)  =  £total(PdHn )  -  fttocailPd  bulk) 
—  —  £totai ( Ho  molecule). 


(6) 


This  definition  yields  £Yorm(PdH)  =  -0.53  e'v  and 
£form(PdHo)  =  +1.34  eV  indicating  a  compound  form¬ 
ing  tendency  of  PdH  and  structural  instability  of  P  Ho, 
again  in  agreement  with  experiment. 

Cohesive  properties  of  the  Pd-H  systems  and  the  ori¬ 
gin  of  the  large  distance  between  hydrogen  atoms  in  bulk 
Pd  can  be  best  understood  by  inspecting  the  electronic 
density  of  states  and  the  corresponding  charge  density. 
In  Fig.  2  we  compare  the  total  density  of  states  of  PdH? 
to  that  of  PdH  (at  the  same  volume)  and  to  bulk  Pd  at 
its  equilibrium  volume.  In  agreement  with  previous  cal¬ 
culations  for  PdH  in  the  bulk  phase4'  and  for  the  surface 
adsorption13  of  II/Pd( 00 1 ) ,  we  find  that  the  binding  be¬ 
tween  H  and  Pd  is  predominantly  covalent,  giving  rise 
to  a  bonding  split-off  state  at  *  —  13  eV  (with  respect 
to  Ep)  in  PdHo  [Fig.  2(b)]  and  at  s  -  8  eV  in  PdH 
[Fig.  2(c)].  To  a  lesser  extent,  hybridization  with  hydro¬ 
gen  also  modifies  the  4d  density  of  states  of  Pd.  With 
increasing  II  concentration,  the  4d  band  is  being  pulled 
below  the  Fermi  level  and  filled,  decreasing  the  contri¬ 
bution  of  d  electrons  in  bonding.  In  the  extreme  case  of 
PdHo,  Fig.  2(c)  suggests  a  transition  to  a  semimetal. 

In  Fig.  3  we  compare  total  charge  densities  of  PdHo, 
PdH,  and  bulk  Pd.  Due  fo  the  large  lattice  constant 
of  Pd,  the  H  atom  and  Ho  molecule  fit  easily  into  the 
Pd  octahedral  cage  [Figs.  3(c)  and  3(d)]  and  fill  the  low 
charge-density  region  inside  [Fig.  3(b)].  Even  at  a  H- 
II  distance  of  0.74  A,  the  separation  betv.ci.-n  II  and  Pd 
atoms  is  close  to  1.7  A.  which  is  a  typical  Pd — II  bond 
length.  Due  to  the  small  size  of  the  H  atom  in  the  octahe¬ 
dral  cage,  the  interaction  potential  with  the  surrounding 
Pd  matrix  is  very  soft  and  gives  rise  to  very  low  II  vibra¬ 
tion  frequencies.  From  the  charge-density  difference  we 
conclude  that  there  is  a  small  net  charge  transfer  from 
Pd  to  II  (a;  0.2  e  near  the  equilibrium  distance  which  is 
typical  for  ail  Pd-H  systems).  From  the  population  uial- 

h 


FIG.  4.  Charge  density  p(r)  corresponding  to  hydrogen-induced  split-oil  states  [at  E  ~  — S  eV  in  Fig.  Jf  b  l  and  E  ~  —  KJ  eV 
in  Fig.  2(c)]  ot  (a)  PdH  and  (b)  PdHj  iuh-u  =  0.74  A).  Both  contour  plots  are  shown  in  the  (LUO)  plane,  which  contains  the 
four  circled  Pd  atoms  in  Fig.  31a).  The  position  of  Pd  atoms  is  marked  by  •:  that  of  II  atoms  by  x.  The  charge  density  is  in 
units  of  1Q— 7  e/a.u.3;  the  ratio  of  two  consecutive  contours  is  p{n  -r  l)/p(n)=1.5. 


ysis  we  conclude  that  this  charge  originates  mainly  from 
the  Pd  5s  band. 

The  inspection  of  the  charge  density  associated  with 
H-induced  split-off  states  in  PdH  and  Pdlls,  shown  in 
Figs.  4(a)  and  4(b),  shows  that  these  states  are  predom¬ 
inantly  of  H  ls-Pd  Ad  character,  and  are  centered  in 
the  middle  of  the  octahedron  cage.  The  large  binding 
energy  of  these  states,  mainly  due  to  the  attractive  H+ 
potential,  increases  with  increasing  hydrogen  concentra¬ 
tion,  which  stabilizes  these  compounds.  This  effect  is  j 
counteracted  by  the  contribution  of  the  H-induced  anti¬ 
bonding  levels.  These  levels,  nearly  empty  at  low  hydro¬ 
gen  concentrations,  are  pulled  to  higher  binding  energies 
and  populated,  which  leads  to  a  destabilization  of  Pd-H 
compounds  at  increasing  hydrogen  concentrations. 

D.  Clean  Pd(001)  and  Pd(110)  surfaces 

Atomic  arrangement  in  the  Pd(001)  and  Pd(UO)  sur-  j 
faces  is  shown  in  Figs.  5(a)  and  5(d).  In  our  calculation,  ( 
these  surfaces  are  modeled  by  slabs  of  finite  thickness. 
Due  to  the  small  value  of  the  Thomas-Fermi  screening 
length  in  metals,  three  to  five  metal  layers  are  sufficient  . 
to  decouple  the  two  sides  of  the  slab.  We  use  a  repeated  ! 
slab  geometry  and  separate  metal  stabs  by  a  “vacuum  j 
slab''  of  sufficient  thickness  (typically  sslO  A)  in  order 
to  avoid  the  surface  charge  density  and  potential  being 
affected  by  neighboring  slabs. 

The  surface  energy  can  be  derived  from  ike  total  en¬ 
ergy  of  the  bulk  and  of  an  n-layer  slab  using 

£?  =  T(fr,olai(slab)  -  n£totaJ(bulk)|.  (7) 

til  our  calculation,  we  used  a  eulotf  ■■nergy  of  t>4  ily  for 
the  Fourier  components  of  the  charge  density  and  po¬ 


tential.  In  order  to  represent  the  charge  decay  in  the 
vacuum  region  accurately,  we  augmented  Pd-based  basis 
functions  (given  in  the  section  on  bulk  Pd)  by  floating 
orbitals  in  the  vacuum  region.  The  good  result  for  the 
Pd(001)  work  function  <t>  =  5.25  eV  (<?  =  5.22  eV  has 
been  observed  in  an  equilibrated  Pd  film'18)  suggests  that 
the  surface  dipole  moment  and  charge  density  have  been 
determined  with  adequate  accuracy.  Our  values  for  the 
surface  energies  of  Pd(001)  and  Pd(110),  based  on  three- 
layer  Pd  slab  calculations,  are  E°(001)  =  0.49  eV/atom 
and  £°(11Q)  =  1.80  eV/atom.  For  computational  rea¬ 
sons,  we  expect  the  true  physical  value  of  the  Pd(110) 
surface  energy  to  be  smaller  than  our  calculated  value 
suggests.49  Nevertheless,  our  LDA  results  are  consistent 
with  an  intuitive  argument  suggesting  that  the  surface 
energy  of  the  more  open  (110)  surface  should  be  larger 
than  that  of  the  close-packed  (00!)  surface. 

Since  no  precise  experimental  data  exist  for  the  surface 
energy  of  these  single-crystal  surfaces,  we  obtain  a  rough 
estimate  using  the  second-moment  approximation  in  the 
tight-binding  technique,50 

E°  a  £^h(bulk)  £[1  -  (Zi/Zb uik)1/2]-  (8) 


Here,  the  summation  is  over  ail  inequivalent  surface 
sites  i,  and  Z,  is  the  corresponding  coordination  num¬ 
ber.  Based  on  this  formula,  and  using  the  experimen-  — 
tal  value  of  the  bulk  cohesive  energy13  £,-oh  —  -5.89  eV, 
the  surface  energies  of  the  (001)  and  (110)  surfaces  of 
Pd.  based  on  Eq.  (8),  are  E)'(OOl)  ~  0.7  eV/atom  and  — 
E°(110)  s  1.1  eV/atom.  Our  LDA  values  are  consistent 

with  the  tight-binding  estimate,  and  also  lie  close  to  es-  _ 

timates  obtained  by  Miedemajl  for  the  surface  energy  ol 
close-packed  surfaces  E°  ~  2.1  J  in':  ~  1.0  eV/atom.  _ 


Differences  in  the  surface  energy  between  these  two 
surfaces  are  reflected  in  the  surface  (or  stab-averaged) 
electronic  density  of  states  (DOS)  Vi E),  shown  in 
Figs.  5(b)  and  5(e).  These  surface  DOS’s  can  be  thought 
of  as  a  weighted  average  of  the  local  DOS  of  the  top¬ 
most  and  second  layer  of  a  Pd  crystal  terminated  with  a 
(001)  or  (110)  surface.  In  agreement  with  the  intuition, 
the  apparent  width  (or  the  second  moment)  of  :V(£)  is 
smaller  for  the  more  open  (110)  surface  than  the  (001) 
surface.  Due  to  the  lower  coordination  of  surface  atoms, 
both  surfaces  show  a  narrower  surface  4 d  band  than  bulk 
Pd  (see  Fig.  2).  Naturally,  the  total  DOS  of  semi-infinite 
samples  equals  the  bulk  DOS  due  to  the  overwhelming 
number  of  bulk  atoms. 

The  total  charge  density  p( r)  of  the  Pd(001)  and  (110) 
surfaces,  shown  in  Figs.  5(c)  and  5(f),  lic°  very  close  to  a 


charge  density  obtained  from  the  superposition  of  atomic 
charge  densities.  We  conclude  that  the  charge  transfer,  - 
the  presence  of  directional  bonding,  and  the  crystal  fieid 
piav  only  a  minor  role  in  the  charge  distrioution. 

IV.  HYDROGEN  ADSORPTION  ON  Pd/001) 

The  adsorption  energy  E^  reflects  the  strength  of  the 
adsorption  bond.  The  value  per  adsorbed  atom  can  be 
obtained  from  the  total-energy  difference  between  the 
adsorption  system,  the  clean  surface,  and  the  free  ad¬ 
sorbate.  For  atomic  hydrogen  adsorbed  on  Pd(001),  we 
obtain 

Eadf©)  =  ^{£totaJ[H/Pd(001)]  -  Etotal[Pd(001)]} 


FIG.  5.  Top  view  of  the  atomic  arram*einent,  surface  (slab-averai;ed)  electronic  density  of  slates  .V(£)  and  charge  density 
p(r)  of  the  PdfOOl)  surface  ((a)-(c)]  and  the  Pd(UO)  surface  ((dMO).  The  contour  plots  of  ,av)  ire  shown  in  the  IU0)  plane 
— —  ( which  contain*  the  (1 10)  direction )  for  the  (001 )  surface  and  in  the  ( Ill)  plane  ( which  contains  the  [  1 1 2]  direction  )  for  the  t  \  1 01 

surface.  The  position  of  Pd  atoms  is  marked  by  •.  The  charge  density  is  in  units  of  10  i  <:/a.u.J.  the  ratio  ol  two  consecutive 
_  contours  is  #>(n  +  l )/p(n 1.5.  6 


FIG.  6.  (a)  Schematic  top  view  of  Pd(001)  and  the  assignment  of  adsorption  sites:  1:  allow  ix  '■  bridge  :  • and  on-top  iC). 

(b)  H  adsorption  energy  £"» d  on  Pd(001i  as  a  function  of  the  adsorption  ncigitt  It.  The  lowest  two  vibration  levels  arc  indicated 
in  the  potential  wells,  (c)  Potential-energy  surface  for  H/Pd(001)  as  a  function  ot  the  hydrogen  position  in  the  surface 
plane,  presented  as  a  3d  plot.  The  position  of  Pd  atoms  is  indicated  by  □  in  (c). 
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(9)y: 


where  the  coverage  @  denotes  the  ratio  between  the  num¬ 
bers  of  adsorbed  and  topmost  layer  substrate  atoms.  In 
our  calculation,  we  consider  hydrogen  to  be  adsorbed  to 
both  sides  of  the  slab  representing  Pd(001).  Then, 


—  Ead(0)  =  2§{  £*««»(H/pdl001)  slabl 

-£'totai[PJ(001)  slab]}  -  £total(II  atom)  . 

_  (10) 

A  schematic  view  of  the  Pd(001)  surface  and  the  high- 
symmetry  adsorption  sites  is  presented  in  Fig.  6(a)  and 

—  the  H-Pd  interaction  potentials  in  these  sites  are  shown  in 
Fig.  6(b).  Table  I  summarizes  the  calculated  equilibrium 
values  of  the  adsorption  bond  energy  £& ,  the  equilibrium 

—  height  (with  respect  to  the  topmost  Pd  layer)  Hq,  and  the 
Harmonic  vibration  frequency  w  for  these  sites. 

A  comparison  of  the  above  adsorption  energy  results, 
—  obtained  for  a  three-layer  Pd  slab,  to  corresponding  five- 
layer  Pd  slab  results,  showed  that  the  differences  lie 
within  the  experimental  error  of  ±0.02  eV.  Similar  agree- 

—  ment  was  obtained  for  the  vibration  frequencies.  Basis 
functions  for  Pd  and  H  have  been  discussed  in  the  pre¬ 
ceding  subsections  on  the  elements.  In  analogy  to  bulk 


Pd.  we  used  an  energy  cutoff  of  64  Ry  in  the  Fourier  ex¬ 
pansion  of  the  charge  density  and  potentials  in  order  to  — 
ensure  complete  convergence  of  the  LDA  spectrum  and 
total  energies.  The  completeness  of  our  basis  set  and  the 
k-point  sampling  of  the  Brilloum  zone  have  been  tested  — 
by  investigating  the  dependence  of  the  total  energy  on 
the  number  of  basis  functions  and  k  points. 

In  order  to  obtain  accurate  values  for  the  heat  of  ad-  — 
sorption  with  a  localized  basis,  we  extended  the  basis  set 
by  including  additional  lloating  orbitals  at  sites  not  oc¬ 
cupied  by  atoms.  The  basis  set  used  for  the  calculation  - 
of  hydrogen  in  hollow  sites  also  included  orbitals  at  all 
bridge  and  on-top  sites.  The  total  energy  for  hydrogen 
in  bridge  or  on-top  sites  has  been  obtained  by  keeping  — 
essentially  the  same  basis  and  “switching'’  the  hydrogen 
ionic  potential  to  the  respective  sites.  For  the  adsorbate- 
free  Pd  slab,  all  basis  sites  at  adsorbate  positions  were  — 
kept  but  carried  zero  potential.  Using  this  procedure, 
we  were  able  to  treat  various  geometries  on  equal  foot¬ 
ing  in  the  orbital  space  and  to  obtain  accurate  energy 
differences  between  the  adsorption  sites. 

Since  our  calculations  are  for  ordered  monolayers,  ad¬ 
sorption  energies  contain  adsorbate-adsorbate  interac-  - 
lions  to  all  orders  and  are  useful  for  the  calculation  of 
adsorption  phase  diagrams.52  As  in  all  LDA  calculations. 


TABLE  I.  Calculated  and  observed  values  of  the  adsorption  bond  energy  adsorption  height 
/to ,  and  vibration  frequency  w  for  H  on  Pd(001). 


Adsorption 

site 

Coverage 

O 

(eV) 

Calc.  Expt. 

Calc. 

ho  (A) 

Expt. 

_  (me  VI 

Calc.  Expt. 

Hollow 

1.0 

—  2.77* 

0.21 

0.30  -  O.U5u 

77  Dl‘ 

Bridge 

2.0 

—2.50 

1.00 

l  Go 

On  -top 

1.0 

-1.36 

1.56 

217 

iWiUUi 


O*  -  hoi 


*  Reference  8. 
b  Reference  12. 
‘Reference  11. 


we  find  our  adsorption  system  slightly  overbound  when 
comparing  to  observed  adsorption  energies.  Previous  cai- 
culations  indicate  t hat  this  overbinding  does  not  signiii- 
cantlv  affect  energy  differences  between  different  adsorp¬ 
tion  sites.  In  Fig.  6(c),  we  shov  a  3d  view  of  the  potential 
energy  £'ad  for  H  along  the  Pdf  001)  surface."3  which  is 
useful  when  discussing  the  surface  diffusion  of  hydrogen. 

When  making  a  comparison  with  the  experiment,  the 
value  of  £ ^  must  also  be  corrected  to  account  for  the 
zero-point-motion  of  the  hydrogen  atom  in  the  potential 
well.  Considering  only  vibrations  perpendicular  to  the 
surface,  and  using  the  calculated  vibration  frequencies  — 
given  in  Table  I,  we  find  the  zero-point-motion  corrected 
values  of  E&i  to  be  —2.88  eV  in  the  hollow  site.  —2.42  eV 
in  the  bridge  site,  and  —1.75  eV  in  the  on-top  site. 

Adsorption  energies  are  related  to  the  measured1*  isos- 
teric  heat  of  adsorption  H ^  and  the  hydrogen  dissocia¬ 
tion  energy  De(Hi)  by 

Hm(H2)  =  -[2£T^(H)-r  D,(H?}]  .  (11) 

Using  our  calculated  value  DS(H:>)  =  4.71  eV,  we  obtain 

=  1-13  eV  for  the  hollow  site.  0.29  eV  for  the  bridge 
site,  and  —0.99  eV  for  the  on-top  site.  Clearly,  hydro¬ 
gen  will  only  adsorb  dissociatively  if  the  net  energy  gain 
Had  >  0.  This  is  the  case  for  the  hollow  and  the  bridge 
sites:  hydrogen  is  not  expected  to  occupy  the  on-top  site. 

For  adsorbate-covered  surfaces,  we  also  detine  a  modi¬ 
fied  surface  energy  (per  surface  Pd  atom)  by 

£,(©)  =  ^{£"totai(H/Pd(001)  slab]  -  n£totaJ( Pd  bulk)} 
-^0  £tot»i(H2  molecule)  .  (12) 

Here,  n  is  the  number  of  Pd  layers  in  a  slab  and  the  total 
energy  of  the  slab  ittotaifH/PdfOOl)  slab]  is  normalized 
to  one  Pd  atom  per  layer.  Using  Eqs.  (7),  (10),  and  ( 11), 
we  can  express  Ea(Q)  as 

£,(©)  =  E°  +  *©[2£«,(©)  +  D«(H2)1 

=  E°  -  l©H*d(H3).  (13) 

In  the  case  if  a  clean  surface,  £,(©  =  0)  corresponds 
to  Eq.  (7).  Our  values  for  the  heats  of  adsorption  sug¬ 
gest  that  hollow-  and  bridge-site  hydrogen  stabilizes  the 
Pd(001)  surface  by  lowering  the  surface  energy.  HydiO- 
gen  in  the  (energetically)  unstable  on-top  site  would  in¬ 
crease  the  surface  energy. 

Our  results  for  the  adsorption  energies,  summarized  in 
Table  I,  establish  that  at  low  coverages,  ©  <  1,  hydrogen 
occupies  preferentially  the  hollow  sites.  This  is  also  sup¬ 
ported  by  the  agreement  between  the  calculated  and  mea¬ 
sured  adsorption  heights13  and  vibration  frequencies^ 

In  the  coverage  range  ©  >  l,  the  additional  hydrogen 
atoms  occupy  new  adsorption  sites.  As  mentioned  above, 
the  on-top  sites  can  be  ruled  out  for  this  purpose.  A 
simultaneous  occupation  of  neighboring  fourfold  hollow 
and  bridge  sites  is  unlikely  because  of  strong  repulsive 
interactions  between  neighboring  adatoms/''*  "8  Tins  is 
also  confirmed  by  our  calculation  for  a  simultaneous  oc¬ 
cupation  of  hollow  and  bridge  sites  on  IM(001),  which 


yields  £\d  =  —1.45  >.*V  at  0  =  3.  corresponding  to  a 
inucii  weaker  bond  than  in  the  on-rop  sue. 

From  the  total-energy  balance  for  ©  >  we  lino  ’hat 
bridge-bonded  hydrogen  at  a  local  coverage  ©  =  .s 
unstable  with  respect  fo  hydrogen  desorption  and  occu-  **• 
pation  ot  noilow  sites.  Inis  energy  balance  is  expected 
to  be  more  favorable  for  the  bridge  sites  in  cases  ot  :ocai  (jhs £ 
coverage  we'l  below  ©  =  2.  We  thus  conclude  that  at 
coverages  0  >  1.  additional  hydrogen  atoms  are  likely  to 
occupy  bridge  sites  at  a  local  coverage  1  <  0  <  2.  or. 
more  likely,  will  diffuse  into  the  bulk.  As  mentioned  in 
Sec.  Ill  C,  the  net  energy  gain  for  the  latter  process  is 
0.53  eV  per  H  atom. 

The  calculated  smail  vaiue  of  the  hydrogen  adsorption 
height  ho  =  0.24  A  given  in  Table  I  for  the  hoilow  site 
is  also  in  very  good  agreement  witn  the  experimental 
result13  of  ho  =  0.30  —  0.05  A  for  the  p|l  x  1)  phase  of  D 
on  Pd.  It  should  be  noted  that  this  adsorption  geometry 
is  very  close  to  the  octahedral  absorption  site  of  H  ;n  bulk 
Pd  (see- Fig.  31.  The  low  value  of  ho  resuits  from  the  iarge 
Pd-H  nearest-neighbor  distance,  which  exceeds  the  sum 
of  the  Pd  and  II  atomic  radii.  This  is  also  the  physical 
origin  of  the  soft  H-Pd  interaction  potential  as  a  function 
of  h  near  the  equilibrium,  as  snow n  in  Fig.  6(b),  which  re¬ 
sults  in  a  very  low  vibration  frequency—’  perpendicular  to 
the  surface.  Our  calculated  vaiue  -  =  7”  meV  agrees  well 
with  the  observed11  frequency  -=t54  meV  for  PdlOOl) 
and  is  also  close  to  the  buik  value  ui  =  68.5  meV  wnicn  has 
been  reported55  for  the  octahedral-site  H  in  bulk  PdH. 

The  H- Pd  interaction  potential  is  very  anharmonic  ana 
shows  a  sharp  increase  for  low  values  of  h  due  to  an  in-  ~~ 
creased  charge-densitv  overlap  between  the  hydrogen  and 
second-layer  Pd  atoms.  Similar  results  for  the  iiollow- 
site  hydrogen  on  Pd(00I)  have  previously  been  obtained 
from  an  effective  medium  caiculationj^y wnich  character-  Pr¬ 
izes  the  adsorption  system  by  a  few  parameters.  From 
our  calculation,  we  find  that  the  Pd  host  electron  density 
at  /•-  varies  by  <  20%  between  the  different  adsorption 
sites.  This  establishes  the  limits  of  validity  of  the  effec¬ 
tive  medium  theory,38  which  averages  out  covalent  effects 
in  cohesion56.  Our  equilibrium  adsorption  geometries  are 
also  consistent  with  a  H — Pd  bond  length  of  i. 6-1.7  A 
in  the  bridge  and  on-top  sites,  which  is  equal  to  the  sum 
of  the  atomic  radii  of  Pd  and  H. 

The  effect  of  potential  anharmomcity  on  the  vibration 
frequency  w  has  been  considered  by  expanding  the  low¬ 
est  vibrational  eigenstates  of  a  general  potential  E*a{h) 
in  terms  of  the  eigenstates  of  a  linear  harmonic  oscillator 
(LHO).  Up  to  nine  LlIC  states  have  been  used  in  the  ex- 
pansion  and  ^  was  determined  from  the  difference  of  the 
two  lowest  vibrational  levels.  The  two  lowest  vibrauonai 
levels  are  shown  in  Fig.  6(b)  together  with  a^i/ii.  In  the 
hollow  site,  this  procedure  yielded  a  vibration  frequency 
that  differed  by  =s8  meV  from  the  harmonic  frequency 
w  =  77  meV.  Our  calculation  of  hydrogen  vibration  fre- 
quencies,  however,  considered  only  II  displacements  per¬ 
pendicular  to  the  surface  at  the  high- symmetry  sites  A 
coupling  to  vibration  modes  parallel  to  the  .surface  is  >-x- 


peeled  to  reuiice  'he  lowest  vibration  frequency  and  bring 
;t  to  closer  agreement  with  the  observed  value-1 

As  seen  in  mg.  bi'b'i  die  hydrogen  adsorption  bond 
strength  decreases  in  rite  sequence  hollow  —  bridge  — 
on- top  site.  In  order  to  understand  this  trend,  we  show 
in  Figs.  7ia  i-7ic  the  surface  slap-averaged)  surface  den¬ 
sity  of  states  Ah  Eh  for  the  different  adsorption  sites.  A 
comparison  with  the  hydrogen-free  Pdf 001)  surface  den¬ 
sity  of  states,  given  ,n  Fig.  5(b),  shows  that  the  most 
prominent  feature  .s  a  hydrogen-induced  sp’it-ofF  state 
below  the  Pd  -id  bam-  in  close  analogy  with  bulk  Pd 
and  PdH  (see  Fig.  0).  The  simple  argument  that  strong 
adsorption  bonus  result  from  large  adsorbate-substrate 
hybridization  m  highly  coordinated  sites  is  supported  by 
comparing  the  energy  difference  ie  between  the  bottom 
of  the  Pd  band  and  the  s, nit-off  H-Pd  bonding  state. 
Tins  quantity  reflects  the  'Hybridization  between  pailr- 
i'um  and  hydrogen.  Inspecting  the  densities  of  states. 


we  lind  no  -niit-olf  states  for  the  on-top  adsorption  site, 
be  ~s  O.J  eV  for  the  bridge  sue.  and  •>  ss  1.1  «*V  rue 
liodow  site.  Of  course,  we  ooserve  a  decrease  >f  y  witn 
increasing  !I-Pd  d’stance.  The  siiailow  potential 
for  hoi  low- site  hydrogen  resuits  from  a  very  siow  variation 
in  tiie  adsorbate  -substrate  hynriuization  a.->  a  function  of 
n.  and  is  retlected  in  relatively  small  changes  of  cy . 

In  Figs.  7(d)-7ff)  we  present  the  totai  charge  den¬ 
sity  of  tiie  [I /Pdf  00 1 1  system  for  tiie  jiiferent  adsorption 
sites  and  h q.  The  changing  environment  of  ’lie  hydro¬ 
gen  atoms  in  tiie  different  adsorption  sites  modifies  not 
oniy  ti.e  depth,  but  also  the  shape  of  the  rl- Pci  interac¬ 
tion  potential,  as  shown  in  F;g.  6(b).  The  charge  den¬ 
sity  for  iioilow-site  hydrogen,  shown  in  F  ig.  7i  d).  exirbits 
a  very  sniail  corrugation,  which  has  been  observed  by 
He  diffraction.’'  With  decreasing  coordination  number. 
Ejdf/i)  resemfc!es  more  a  pairwise  interaction  potential 
that  gives  rise  to  an  increasing  vibration  frequency  _ 


FIG.  7  iurlace  i  slab-averaged)  electronic  density  of  slates  VI /S' l  lor  l!,Pdi"Ul)  ,n  tin-  hollow  ,ri.  imugc-  lie  urn  on-’on 
-lie  lc).  Charge  density  ,n  r  I  for  11/  PdfOU  l )  in  the  hollow  id),  i  ridge  mi.  .ml  on- 1  op  - 1  ■  -  l  l ,  -now  n  n  i  In-  ■  l  l  n  l  ;n.i  ne  n  imn 
oniailis  l  lie  Alt'1  direction).  The  position  of  Pd  ilums  is  marked  Uv  •  mo  ili.u  ->  l  !(  uom-  <  ihe  ih.uge  Ucumh  -  m 
units  of  1 1)“ J  •■/  .t.ii. "  the  ratio  of  two  consecutive  i  on  tours  is  ,u  "  —  1  i/Vi  «  .  =  i  ‘ 


FIG.  8.  (a)  Difference  density  of  c'ates  Ai VIE),  (b)  difference  charge  density  Aphu  and  fcj  charge  densiiv  .associated  with 

die  H-induced  split-off  state  [at  E  ~  —7  eV  in  Fig.  7(a)]  fo"  II/Pdf001)  in  the  lioilc  .  site.  The  position  o’  Pd  atoms  is  m.arKed 

_  by  •  and  that  of  H  atoms  by  x.  The  charge  density  is  in  units  of  10_i  e/a..  3.  In  (b),  tlie  consecutive  contours  are  separated 

by  ID-*  e/a.u.3 ,  and  in  (cj  the  ratio  of  two  consecutive  contours  is  p\n  —  1  )j p\  ,j)=1.5. 


The  work  function  depends  very  sensitively  on  the  sur¬ 
face  charge  distribution,  winch  establishes  the  surface 

_ dipole.  Our  calculated  work  function  for  the  hollow  site 

(0  =  1)  o  =  5.53  eV  has  increased  by  0.3  eV  with  respect 
to  the  hydrogen-free  surface;  a  similar  increase  of  O.il  eV 
_  has  been  determined  experimentally  for  this  surface.3 

To  establish  the  nature  of  the  hydrogen-palladium  sur¬ 
face  bond,  we  determine  the  difference  density  of  states 

-  A,V( E)  =  .VH/p<hooi)(£)  ~  -Vpdiooi )f£)  and  the  charge- 

density  difference  Ap  £:  pH/Pdiooi)  -  Ppdioon  -  PH  for 
the  hollow-site  hydrogen.  An  inspection  of  the  difference 

_  <’  insity  of  states  for  the  hollow  site,  given  in  Fig.  S(a), 

shows  that  the  main  effect  of  hydrogen  is  to  induce  a 
split-off  state  below  the  bulk  Pd  band,  and — to  a  lesser 

_  extent — modify  the  Pd  Ad  band.  An  analysis  of  the 

charge-der  ty  diffe-ence,  given  in  Fig,  8(b),  shows  a 
charge  transfer  of=:0.2  electrons  towards  hydrogen,  stem- 

_  mtng  mainly  from  Pd  atoms  in  the  topmos'  layer.  The 

nature  of  the  II- Pd  bonding  mechanism  can  be  seen  from 
the  charge  density  associated  with  the  H-Pd  split-off 


state  at  s;  I  eV  below  the  Fermi  level  in  Fig.  8(a),  shown 
in  Fig.  8(c).  The  covalent  bonds  result  main’y  from  the 
hybridization  of  II  Is  and  Pd  Ad  states  and  are  mostly 
localized  in  the  topmost  surface  layer.  These  states  are 
very  close  in  character  to  the  11-induced  spiit-off  states 
in  bulk  PdH  and  »  dlls,  which  are  analyzed  in  Fig.  4. 

V.  HYDROGEN  ADSORPTION  ON  Pd(110) 

In  this  section,  we  will  discuss  hydrogen  adsorption  on 
Pd(ilO)  along  the  lines  of  H/?d(00I).  When  presenting 

numerical  results,  we  will  use  obvious  modifications  of  _ 

the  expressions  for  the  adsorption  bond  energy  Ead  and 
surface  energy  E,(0)  given  in  Eqs.  (9),  (10)  and  (12). 

As  illustrated  in  Figs.  3(d)  and  9(a),  Pd(lIO)  is  a  very  _ 

open  surface.  This  is  the  origin  of  the  large  surface  energy 
value  discussed  in  Sec.  Ill  D,  but  also  raises  the  question 

of  surface  relaxations.  Figure  9(a)  gives  a  side  view  of _ 

a  relaxed  Pd(I10)  surface.  The  vertical  relaxation,  indi¬ 
cated  in  Fig.  9(a),  is  defined  by  A =  (dl2/d?3  —  1), 


FTG.  'J.  <a|  Schematic  side  vie'"  of  the  Pdf  110)  surface  (b)  Surface  energy  as  a  (unction  oi  the  vertical  relaxation  Akij  of 
the  lirsl  inlerlaver  spacing  in  Pdf  l  10).  The  solid  line  corresponds  to  a  clean  Pdf  1  10)  surface  and  the  dashed  line  to  a  II-<  ovned 
|0=l,  hollow  sitel  surface.  j 


where  d?,  =  1.38  A  is  the  bulk  interlayer  spacing.  Our  re¬ 
sults  for  the  surface  energy  of  clean  Pd(  110),  given  by  the 
solid  line  in  Fig.  9(b),  indicate  a  contraction  d\2  =  —4% 
of  the  topmost  layer  spacing.  In  view  of  the  fact  that 
we  did  not  consider  multilayer  relaxations,  this  value 
agrees  well  with  the  LEED  result14  d12  =  —5.1%  and 
with  the  general  trend  observed  and  calculated  for  sur¬ 
faces  of  late  transition  metals. l4,5‘  It  is  conceivable  that 
a  separate  consideration  of  multilayer  relaxations  could 
enhance  Adi2  (by  reducing  the  repulsion  between  the 
topmost  and  the  third  layer).58 

In  anaiogy  to  H/Pd(001),  H  adsorption  lowers  the  sur¬ 
face  energy  by  saturating  Pd  dangling  bonds  and  con¬ 
sequently  reverses  the  contraction.  This  is  shown  by 
the  dashed  line  in  Fig.  9(b),  which  corresponds  to  a 
hydrogen-covered  surface.  In  our  calculation,  we  as¬ 
sumed  the  H  atoms  to  be  adsorbed  in  the  hollow  site 
(see  Fig.  10)  and  kept  the  distance  between  hydrogen  and 
second-layer  Pd  atoms  fixed.  For  a  monolayer  coverage 
of  hydrogen,  our  results  indicate  a  6%  surface  expansion. 
This  expansion  is  reminiscent  of  the  value  observed  in 
bulk  Pd  upon  hydrogen  uptake  and  reflects  similarities 
in  the  Pd — H  bond. 

As  shown  in  Fig.  9(b),  the  Pd-Pd  interplane  interaction 
potential  has  a  large  anharmonicity,  which  has  also  been 
observed  on  the  chemically  similar  Ni(  110)  surface.39  The 
interplane  force  constant  c  =  12.3  eV/A2  at  the  Pd(110) 
surface  is  somewhat  smaller  than  the  bulk  force  constant 
c  =  15.4  eV /A2  corresponding  to  a  compression  along  the 
[110]  direction.  Hydrogen  adsorption  on  Pd(110)  lowers 
the  interplane  force  constant  to  c  =  10.3  eV/A2.  A  sim¬ 
ilar  hydrogen-induced  softening  can  be  observed  in  the 
bulk  (see  Fig.  1),  which  corresponds  to  a  decrease  of  the 
Pd  bulk  modulus  in  the  hydride  phase. 

A  schematic  top  view  of  the  Pd(llO)  surface  and  the 
high-svmmetry  adsorption  sites  is  presented  in  Fig.  10(a) 


and  the  H-Pd  interaction  potentials  in  these  sites  are 
shown  in  Fig.  10(b).  All  calculations  have  been  per¬ 
formed  in  an  analogous  way  to  the  Pdf 001)  surface  de¬ 
scribed  in  the  preceding  section,  for  five-layer  slabs  con¬ 
sisting  of  three  palladium  layers  and  two  hydrogen  lay¬ 
ers.  Table  II  summarizes  the  calculated  equilibrium  val¬ 
ues  of  the  adsorption  bond  energy  ald,  the  equilibrium 
height  (with  respect  to  the  topmost  Pd  layer)  ho,  and 
the  harmonic  vibration  frequency  w  for  these  sites.  A  3d 
view  of  the  potential  energy  £* d  for  H  along  the  Pd(110) 
surface53  is  shown  in  Fig.  10(c). 

When  making  comparison  with  the  experiment,  the 
value  of  Ead  has  also  to  be  corrected  to  account  for  the 
zero-point-motion  of  the  hydrogen  atom  in  the  potential 
well.  Considering  only  vibrations  perpendicular  to  the 
surface,  and  using  the  calculated  vibration  frequencies  ^ 
given  in  Table  II,  we  find  the  zero-point-motion  corrected 
values  of  E^d  t0  he  —3.38  eV  in  the  long-bridge  site, 
—  2.98  eV  in  the  hollow  site,  —3.00  eV  in  the  short-bridge 
site,  and  —2.17  eV  in  the  on-top  site.  Due  to  numeri¬ 
cal  difficulties  similar  to  the  calculation  of  the  Pd(110) 
surface  energy,  these  calculated  values  of  £ad  probably 
exaggerate  the  strength  of  the  adsorption  bond.49 

Using  our  calculated  value  D ,(//2)  =  4.71  eV,  we  ob¬ 
tain  =  2.11  eV  in  the  long-bridge  site.  1.43  eV  in  the 
hollow  site,  1.39  eV  in  the  short-bridge  site,  and  —0.11  eV 
in  the  on-top  site.  Cle.  hydrogen  will  only  adsorb  dis- 
sociatively  if  the  net  energy  gain  Had  >  0.  This  is  the 
case  for  the  long-bridge,  hollow,  and  short-bridge  sites. 
Similar  to  the  Pd(001)  surface,  hydrogen  is  not  expected 
to  occupy  the  on-top  site. 

Using  Eq.  (13)  and  these  values  for  the  heats  of  adsorp¬ 
tion.  we  can  determine  the  surface  energy  of  hydrogen- 
covered  Pd(110).  Our  results  suggest  that  hydrogen  in 
long-bridge,  hollow,  and  short-bridge  sites  stabilizes  the 
Pd(110)  surface  by  lowering  the  surface  energy.  Hydro- 


FIG.  10.  (al  Schematic  top  view  of  Pdf  1 10)  and  tlte  assignment  of  adsorption  sues:  hollow  |xt.  !oug-iir»lge  i»t.  short- 
bridge  lo),  and  on-top  (C'.  i  b )  il  adsorption  energy  on  Pdf  1 101  as  a  function  ol  the  adsorption  height  h.  The  lowest 
two  vibration  levels  are  indicated  in  the  potential  wells,  (c)  Potential-energy  surface  E.o  for  U/P  lflHM  as  a  function  of  the 
hydrogen  position  in  the  surface  plane,  presented  .u»  a  J d  plot.  Tiie  position  ol  Pd  atoms  is  indicated  hv  Q  ui  icl.  >. 


TABLE  II.  Calculated  values  of  the  adsorption  bond  energy  £ia ,  adsorption  height 

/<o.  .uul 

vibration  frequency  ui  foi 

r  H  on  Pdf  1 10). 

Adsorption 

Coverage 

£«•  leV) 

ho  (A) 

-  I  nic-V) 

site 

© 

Calc. 

Caic. 

Calc. 

X  W»u.vi 

Long-bridge 

1.0 

-.1.41 

0.05 

O  i 

Hollow 

1.0 

-3.07 

0.35 

159 

Short-b  ridge 

1.0 

-3.05 

1.0 

100 

On-top 

1.0 

-2.30 

1.53 

253 

J 

OM.— fpp 

..  ’ '  ' 

gen  in  the  (energetically)  unstable 

on-top  site  would 

in-  hydrogen 

occupies  preferentially 

the  trough 

sites,  in 

crease  the  surface  energy. 

agreement 

with  experimental  data  obtained 

by  LEED 

Our  results  for  the  adsorption  energies,  summarized  and  LEIS. 14,15  Similar  preferential  adsorption  sites  have 

in  Table  II,  establish  that  at  low  coverages.  0  <  1,  previously  been  found  using  the  embedded-atom"6  and 


FIG.  11.  Surface  (slab-averaged)  electronic  density  of  states  N(E)  for  !I/Pd(110)  in  (a)  the  hollow,  (b)  long-bridge,  (c) 
short-bridge,  and  (d)  on-top  site.  Charge  density  p( r )  for  H/Pd(UO)  in  (e)  the  hollow,  (f)  long-bridge,  (g)  shon-bridgc.  .uul 
(h)  on-top  site.  The  contour  plots  (e)  and  (li)  are  shown  in  the  (111)  plane  (which  contains  the  (112]  direction  I  and  the  plots  if) 
and  (g)  are  given  in  the  (001)  plane  (which  contains  the  surface  (110|  direction  |.  The  position  of  Pd  atoms  is  marked  bv  •  and 
Uiai  of  II  atoms  by  x .  The  charge  density  is  in  units  of  IQ-2  e/a.u.2 ,  the  ratio  of  two  consecutive  contours  is  ,n  n  1  )/,u  »  |=  i  .2. 


•.‘mbedded-clusterjr  method.  At  higher  coverages.  0  >  i. 
new  adsorption  sues  are  occupied.  Eventually,  at  large 
coverages  exceeding  0  ~  1.5,  a  transition  to  a  streaky 
I  1  x  21  phase  involving  substrate  reconstruction  has  been 
oDserved  21--- 

Basea  on  a  careful  analysis  of  LEED  (Ret.  14)  and 
He  diffraction10  data,  the  quasithreefold  site  located  in 
troughs  between  the  hollow  and  the  short-bridge  sites 
shown  in  Fig.  10(ai]  has  been  suggested  as  the  preferen¬ 
tial  adsorption  site  for  H  in  the  (2  x  1)  phase  observed  at 
low  coverages  0  <  1.  Due  to  similarities  in  the  chemical 
environment,  we  expect  the  nature  of  the  H — Pd  bond 
in  the  quasithreefold  site  to  be  very  similar  to  the  more 
symmetric  bridge  and  hollow  sites.  We  have  performed 
no  calculations  for  the  quasithreefold  sue  since  the  com¬ 
putational  requirements  for  this  low-symmetry  sue  are 
very  high. 

Our  equilibrium  adsorption  geometries  are  consistent 
with  a  H — Pd  bond  length  of  2.0  A  in  the  twofold  co¬ 
ordinated  long-bridge  and  short-bridge  sues.  This  is  in 
excellent  agreement  with  LEED  results*-1  for  the  H  —  Pd 
nearest-neighbor  distance  of  2.0±0.1  A  for  H/Pd(110). 
For  the  singly  coordinated  on-top  and  hollow  sues,  the 
H — Pd  bond  length  ranges  between  1.5  and  1.7  A.  which 
is  consistent  with  the  sum  of  atomic  radii  of  H  and  Pd. 

At  low  hydrogen  coverages  0  <  1.  two  vibration 
modes  have  been  observed1, •l£t  on  Pd(110)  single  crys¬ 
tals  at  rslOO  and  a  120  meV.  Assuming  adsorption  in  the 
quasithreefold  site  of  Pd(llQ),  based  an  analogy  to 
similar  modes  observed60  for  II /Pdf  111),  the  higher- 
frequency  mode  has  been  assigned  to  a  vibration  per¬ 
pendicular  to  the  plane  of  the  three  Pd  neighbors  of  H, 
and  the  lower  mode  to  a  motion  parallel  to  this  plane. 
While  the  calculation  of  these  vcry-low-symmetry  sites 
and  the  normal  modes  of  II/Pd(  1 10)  is  computationally 
too  demanding,  our  values  for  u  in  the  chemically  similar 
hollow  and  short-bridge  sites,  presented  in  Table  II,  lie 
in  the  same  range  as  the  experimental  values.  Our  calcu¬ 
lated  value  w=57  meV  tor  the  long-bridge  site  seems  to 
match  perfectly  with  a  low-frequency  mode  -=53  meV, 
which  has  been  observed  on  Pd  powder  and  assigned  to 
subsurface  sites.61 

As  seen  in  Fig.  10(b),  the  hydrogen  adsorption  bond 
strength  decreases  in  the  sequence  long-bridge  —  hoilow 


—  short-bridge  —  on-top  site.  In  order  to  understand 
this  trend,  we  show  m  Figs.  1 1(a)— 1 1(c)  the  surtace  i  slab- 
averaged)  surface  density  of  states  ,V(  E)  for  the  differ¬ 
ent  adsorption  sites.  A  comparison  with  the  hydrogen- 
free  Pdf  110)  surface  density  of  states,  given  in  Fig  5(e), 
shows  that  the  most  prominent  feature  is  a  hydrogen- 
induced  split-off  state  below  the  Pd  4 d  band,  in  close 
analogy  with  the  hydrogen-covered  Pd(001)  surface  [see 
Figs.  7(a)-7(d)  ].  In  analogy  with  the  Pd(001)  results, 
we  can  correlate  strong  adsorption  bonds  with  a  large 
adsorbate-substrate  hybridization  in  highly  coordinated 
sites.  This  is  supported  by  comparing  the  binding  ener¬ 
gies  of  the  split-off  state  with  respect  to  the  bo  i  of  the 
Pd  band,  which  reflects  the  hybridization  betw_.n  palla¬ 
dium  and  hydrogen.  Inspecting  the  densities  of  states,  ive 
find  the  binding  energy  of  the  split-off  state  to  be  largest 
in  the  long-bridge  site  and  smallest  in  the  on-top  site. 

In  Figs.  1 1(e)—  1 1  f  h )  we  present  the  total  charge  den¬ 
sity  of  the  H/Pd(  110)  system  for  the  different  adsorption 
sites  for  the  equilibrium  height  h0.  The  changing  environ¬ 
ment  of  the  hydrogen  atoms  in  the  different  adsorption 
sites  modifies  not  only  the  depth,  but  also  the  shape  of 
the  H-Pd  interaction  potential,  as  shown  in  Fig.  10(b). 
As  illustrated  in  Fig.  11(e)  for  the  case  of  a  hollow  site, 
hydrogen  prefers  to  be  “buried”  in  the  trough  sites  caus¬ 
ing  a  smoothening  of  the  surface  charge  density  and  a  low 
vibration  frequency.  With  decreasing  coordination  num¬ 
ber.  Erui{h)  resembles  more  a  pairwise  interaction  poten¬ 
tial  that  gives  rise  to  an  increasing  vibracion  frequency  ^ 
[see  the  two  lowest  levels  in  Fig.  10(b)], 

To  establish  the  nature  of  the  hydrogen-palladium  sur¬ 
face  bond,  we  determine  the  difference  density  of  states 
A(V(£)  =  .VH/pd(Ut»(£)  -  iVPd(110)(E)  and  the  charge- 
density  difference  Ap  =  PH/Pddio)  -  pPdiuot  -  Pu  for 
the  hollow-site  hydrogen.  An  inspection  of  the  difference 
density  of  states  for  the  hollow  site,  given  in  Fig.  12(a), 
shows  tiiat  the  main  effect  of  hydrogen  is  to  induce  a 
split-off  state  below  the  buik  Pd  band,  and — to  a  lesser 
extent — modify  the  Pd  4 d  band.  An  analysis  of  the 
charge-density  difference,  given  in  Fig.  12(b),  shows  a 
charge  transfer  of  ss0.2  electrons  towards  hydrogen,  stem¬ 
ming  mainly  from  Pd  atoms  in  the  underlying  layer. 
The  nature  of  the  H-Pd  bonding  mechanism  can  be  seen 
from  the  charge  density  associated  with  the  H-Pd  split- 


FIG.  12.  (a)  Difference  density  of  states  A.V(£),  (b)  difference  charge  dvusitv  A/hi'l.  and  icl  ■  haigc  deiiMi 

with  t  lie  U-induced  split-off  slate  (at  E  s:  -  T  eV  in  Fit,.  11(a)]  for  II  /  Pd  I  1 10 )  in  the  hollow  site.  The  pomiioh  .u 
marked  by  •  and  that  of  ![  atoms  by  x.  ".  ue  charge  density  is  in  units  nl  1U“‘  e/a.ii."  In  ibi.  i  lie  coiixcinivc 
separated  by  10”1  e/a.u.1,  and  in  (c)  the  ratio  ol  two  consecutive  contours  is  ,i\  n  -+•  l  )//U  '*  l=  1 
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off  state  at  as  7  eV  below  the  Fermi  level  in  Fig.  12(a), 
shown  in  Fig.  12(c).  The  covalent  bonds  result  mainly 
from  the  hybridization  of  H  Is  and  Pd  4d  states  of  second- 
layer  substrate  atoms  neighboring  the  adsorbate.  These 
states  are  very  close  in  character  to  the  H-induced  split- 
off  states  in  H/Pd(001),  which  are  analyzed  in  Fig.  8. 

VI.  SUMMARY  AND  CONCLUSIONS 

In  summary,  we  used  the  ab  initio  pseudopotential  lo¬ 
cal  orbital  method34  to  determine  the  interaction  of  hy¬ 
drogen  with  the  (001)  and  (110)  surfaces  of  Pd  as  well  as 
bulk  Pd.  Our  calculations  for  the  the  H-Pd  system  have 
been  completed  by  calculations  for  the  H?  molecule,  the 
bulk,  and  the  clean  (001),  and  (110)  surfaces  of  Pd. 

We  found  that  the  preferential  adsorption  sites  for  hy¬ 
drogen  are  the  trough  sites  on  Pd(110)  and  the  fourfold  : 
hollow  site  on  Pd(001)  which  corresponds  to  the  octahe¬ 
dral  site  in  bulk  Pd.  The  formation  of  a  H — Pd  bond 
is  accompanied  by  a  transfer  of  xO.2  electrons  from  Pd 
to  hydrogen.  The  signature  of  the  H — Pd  bond,  which 
has  a  covalent  H  ls-Pd  Ad  character  and  involves  only 
the  nearest  neighbors  of  the  hydrogen  atom,  is  a  split-off 
state  located  a  few  eV  below  the  lower  Pd  Ad  band  edge. 
The  presence  of  hydrogen  in  the  bulk  and  at  the  surface  of 
Pd  causes  a  shift  of  the  Pd  Ad  band  to  higher  binding  en¬ 
ergies.  For  a  nearly  full  Ad  band,  this  increased  band  fill¬ 
ing  results  in  decohesion  and  softer  force  constants  both 


at  the  surface  and  in  the  bulk.  The  calculated  and  ob¬ 
served  lattice  constant  increase  of  s;6%  in  bulk  Pd  upon 
hydrogen  uptake  in  the  high-density  phase  corresponds  to 
a  small  hydrogen-induced  expansion  of  otherwise  slightly 
contracted  Pd  surfaces.  In  line  with  the  general  trend  ob¬ 
served  in  chemisorption  systems,  lower-coordinated  ad¬ 
sorption  sites  are  energetically  less  favorable  for  hydro¬ 
gen.  With  decreasing  coordination  number,  the  H — Pd 
bond  length  decreases  from  «2.Q  to  ssl.5  A  and  the  H-Pd 
interaction  potential  becomes  narrower,  thus  resulting  in 
higher  vibration  frequencies. 
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